Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics/Quantum Chemistry and Polaron Model Study

The dynamics of pigment-pigment and pigment-protein interactions in light-harvesting complexes is studied with a novel approach which combines molecular dynamics (MD) simulations with quantum chemistry (QC) calculations. The MD simulations of an LH-II complex, solvated and embedded in a lipid bilayer at physiological conditions (with total system size of 87,055 atoms) revealed a pathway of a water molecule into the B800 binding site, as well as increased dimerization within the B850 BChl ring, as compared to the dimerization found for the crystal structure. The fluctuations of pigment (B850 BChl) excitation energies, as a function of time, were determined via ab initio QC calculations based on the geometries that emerged from the MD simulations. From the results of these calculations we constructed a time-dependent Hamiltonian of the B850 exciton system from which we determined the linear absorption spectrum. Finally, a polaron model is introduced to describe quantum mechanically both the excitonic and vibrational (phonon) degrees of freedom. The exciton-phonon coupling that enters into the polaron model, and the corresponding phonon spectral function are derived from the MD/QC simulations. It is demonstrated that, in the framework of the polaron model, the absorption spectrum of the B850 excitons can be calculated from the autocorrelation function of the excitation energies of individual BChls, which is readily available from the combined MD/QC simulations. The obtained result is in good agreement with the experimentally measured absorption spectrum.Comment: REVTeX3.1, 23 pages, 13 (EPS) figures included. A high quality PDF file of the paper is available at http://www.ks.uiuc.edu/Publications/Papers/PDF/DAMJ2001/DAMJ2001.pd

A Comparative Molecular Dynamics Study of Methylation State Specificity of JMJD2A

Histone modifications have great importance in epigenetic regulation. JMJD2A is a histone demethylase which is selective for di- and trimethyl forms of residues Lys9 and Lys36 of Histone 3 tail (H3K9 and H3K36). We present a molecular dynamics simulations of mono-, di- and trimethylated histone tails in complex with JMJD2A catalytic domain to gain insight into how JMJD2A discriminates between the methylation states of H3K9. The methyl groups are located at specific distances and orientations with respect to Fe(II) in methylammonium binding pocket. For the trimethyllysine the mechanism which provides the effectual orientation of methyl groups is the symmetry, whereas for the dimethyllysine case the determining factors are the interactions between methyllysine head and its environment and subsequently the restriction on angular motion. The occurrence frequency of methyl groups in a certain proximity of Fe(II) comes out as the explanation of the enzyme activity difference on di- and tri-methylated peptides. Energy analysis suggests that recognition is mostly driven by van der Waals and followed by Coulombic interactions in the enzyme-substrate interface. The number (mono, di or tri) and orientations of methyl groups and water molecules significantly affect the extent of van der Waals interaction strengths. Hydrogen bonding analysis suggests that the interaction between JMJD2A and its substrates mainly comes from main chain-side chain interactions. Binding free energy analysis points out Arg8 as an important residue forming an intra-substrate hydrogen bond with tri and dimethylated Lys9 of the H3 chain. Our study provides new insights into how JMJD2A discriminates between its substrates from both a structural and dynamical point of view

Molecular Recognition of H3/H4 Histone Tails by the Tudor Domains of JMJD2A: A Comparative Molecular Dynamics Simulations Study

Background: Histone demethylase, JMJD2A, specifically recognizes and binds to methylated lysine residues at histone H3 and H4 tails (especially trimethylated H3K4 (H3K4me3), trimethylated H3K9 (H3K9me3) and di, trimethylated H4K20 (H4K20me2, H4K20me3)) via its tandem tudor domains. Crystal structures of JMJD2A-tudor binding to H3K4me3 and H4K20me3 peptides are available whereas the others are not. Complete picture of the recognition of the four histone peptides by the tandem tudor domains yet remains to be clarified. Methodology/Principal Findings: We report a detailed molecular dynamics simulation and binding energy analysis of the recognition of JMJD2A-tudor with four different histone tails. 25 ns fully unrestrained molecular dynamics simulations are carried out for each of the bound and free structures. We investigate the important hydrogen bonds and electrostatic interactions between the tudor domains and the peptide molecules and identify the critical residues that stabilize the complexes. Our binding free energy calculations show that H4K20me2 and H3K9me3 peptides have the highest and lowest affinity to JMJD2A-tudor, respectively. We also show that H4K20me2 peptide adopts the same binding mode with H4K20me3 peptide, and H3K9me3 peptide adopts the same binding mode with H3K4me3 peptide. Decomposition of the enthalpic and the entropic contributions to the binding free energies indicate that the recognition of the histone peptides is mainly driven by favourable van der Waals interactions. Residue decomposition of the binding free energies with backbone and side chain contributions as well as their energetic constituents identify the hotspots in the binding interface of the structures. Conclusion: Energetic investigations of the four complexes suggest that many of the residues involved in the interactions are common. However, we found two receptor residues that were related to selective binding of the H3 and H4 ligands. Modifications or mutations on one of these residues can selectively alter the recognition of the H3 tails or the H4 tails

Direct observation of topoisomerase IA gate dynamics

Type IA topoisomerases cleave single-stranded DNA and relieve negative supercoils in discrete steps corresponding to the passage of the intact DNA strand through the cleaved strand. Although type IA topoisomerases are assumed to accomplish this strand passage via a protein-mediated DNA gate, opening of this gate has never been observed. We developed a single-molecule assay to directly measure gate opening of the Escherichia coli type IA topoisomerases I and III. We found that after cleavage of single-stranded DNA, the protein gate opens by as much as 6.6 nm and can close against forces in excess of 16 pN. Key differences in the cleavage, ligation, and gate dynamics of these two enzymes provide insights into their different cellular functions. The single-molecule results are broadly consistent with conformational changes obtained from molecular dynamics simulations. These results allowed us to develop a mechanistic model of interactions between type IA topoisomerases and single-stranded DNA

Mechanisms of Hybrid Oligomer Formation in the Pathogenesis of Combined Alzheimer's and Parkinson's Diseases

Background: Misfolding and pathological aggregation of neuronal proteins has been proposed to play a critical role in the pathogenesis of neurodegenerative disorders. Alzheimer’s disease (AD) and Parkinson’s disease (PD) are frequent neurodegenerative diseases of the aging population. While progressive accumulation of amyloid b protein (Ab) oligomers has been identified as one of the central toxic events in AD, accumulation of a-synuclein (a-syn) resulting in the formation of oligomers and protofibrils has been linked to PD and Lewy body Disease (LBD). We have recently shown that Ab promotes a-syn aggregation and toxic conversion in vivo, suggesting that abnormal interactions between misfolded proteins might contribute to disease pathogenesis. However the molecular characteristics and consequences of these interactions are not completely clear. Methodology/Principal Findings: In order to understand the molecular mechanisms involved in potential Ab/a-syn interactions, immunoblot, molecular modeling, and in vitro studies with a-syn and Ab were performed. We showed in vivo in the brains of patients with AD/PD and in transgenic mice, Ab and a-synuclein co-immunoprecipitate and form complexes. Molecular modeling and simulations showed that Ab binds a-syn monomers, homodimers, and trimers, forming hybrid ringlike pentamers. Interactions occurred between the N-terminus of Ab and the N-terminus and C-terminus of a-syn. Interacting a-syn and Ab dimers that dock on the membrane incorporated additional a-syn molecules, leading to th

A Distributed Dynamic Load Balancer for Iterative Applications

For many applications, computation load varies over time. Such applications require dynamic load balancing to improve the performance. Centralized load balancing schemes, which collect global information and compute the decisions at a central location, are not scalable. In contrast, fully distributed strategies are scalable but typically do not produce a balanced work distribution as they tend to consider only local information. This paper describes a fully distributed algorithm for load balancing that uses partial information about the global state of the system to make good load balancing decisions. This algorithm, referred to as Grapevine, consists of two stages: global information propagation using a lightweight algorithm inspired by epidemic [26] algorithms, and work unit transfer using a randomized algorithm. The proposed algorithm is scalable and can be tuned to optimize for either cost or performance. We substantiate our claims by providing analysis of the algorithm, detailed simulation and performance comparison on adaptive mesh refinement and molecular dynamics applications. We demonstrate the effectiveness of Grapevine for adaptive mesh refinement and molecular dynamics on up to 131,072 cores of BlueGene/Q. 1

Agents: an Undistorted Representation of Problem Structure

It has been observed that data-parallel languages are only suited to problems with "regular" structures. This observation prompts a question: to what extent are other parallel programming languages specialized to specific problem structures, and are there any truly general-purpose parallel programming languages, suited to all problem structures? In this paper, we define our concept of "problem structure". Given this definition, we describe what it means for a language construct to "directly reflect" a problem structure, and we argue the importance of using a language construct which reflects the problem structure. We describe the difficulties that arise when the language construct and the problem structure do not fit each other. We consider existing language constructs to identify the structures they fit, and we note that language constructs are often designed with little regard for such generality. Finally, we describe a parallel language construct which is designed specifically wit..